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The paper investigates the use of active control for stability 
augmentation of passive gravity gradient satellites. The reaction jet 
method of control is the main interest. Satellite nonrigidity is emphasized. 
The reduction in the Hamiltonian H is used as a control criteria. The 
velocities , relative to local vertical, of the jets along their force 

axes are shown to be of fundamental significance. A basic control scheme 
which satisfies the H reduction criteria is developed. Each jet ©< is fired 
when its nj ~uc< becomes appropriately large. The jet is de-energized when 
't’T reaches zero. Firing constraints to preclude orbit alteration may be 
needed. Control is continued until H has been minimized. This control 
policy is investigated using impulse and rectangular pulse models of the 
jet outputs. The impulse model leads to a simple equation for the optimal 
instantaneous control magnitude. This model, however, is difficult to employ 
in nonrigid satellite applications due to basic problems that are discussed. 
The study using the rectangular pulse model includes (1) development of 
general equations for the pulse duration which yields 0 and (2) investi- 

gation of a strategy to override this criteria when necessary to prevent 
large structural vibrations. 
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NOMENCLATURE 

semi-major axis of orbit 

orbit eccentricity 
unit vectors along axes of frame b 
unit vector parallel to orbit pole 
unit vector parallel to local vertical 
unit vector along force axis of jet <=K 
torque vector 

inertia dyadic about satellite's c.m. 
arbitrary constant 
satellite mass 

component of N due to viscous damping 
distance from planet's c.m. to satellite's c.m. 
kinetic energy component of zeroth order in 4 
time at which a control operation is begun 
dynamic potential 
potential energy 

the component of V due to internal stiffness 
vector from satellite's c.m. to jet o<r 
true anomaly angle 
gravitational constant 

angular velocity vector of frame b relative to frames R and r 
true anomaly angular velocity vector 
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INTRODUCTION 


This paper originated in a study performed for the RAE-B satellite. 
RAE-B will have four 750' gravity gradient booms, a passive damper system, 
and six freon reaction jets. Ihe jets will be used for spin vector control 
early in the mission before deployment of the booms and capture by the 
gravity field. The purpose of the study was to investigate the feasibility 
of a suggestion, raised by GSFC personnel, to also use the jets for attitude 
stability augmentation during or after boom deployment if an emergency situ- 
ation is encountered or if unsatisfactory performance of the damper system 
is experienced. Feasibility hinged largely on whether or not properly 
timed force pulses could alleviate attitude librations or tumbling on RAE-B 
without generating intolerably large boom vibrations. The difficulty arises 
because the long booms will be extremely flexible, and the force outputs of 
the jets will be comparatively large. 

The material which has been included in the paper is not limited to 
satellites of the RAE-B configuration. Instead, it is directly applicable 
to virt ually any satellite which is stabilized by the passive gravity 
gradient method and includes one or more reaction jets. While the analytical 
methods which are employed do have some potentiality for extension to other 
types of control, the study is devoted almost entirely to techniques, 
particularly reaction jets, which provide damping or control by the appli- 
cation of external forces to the system. 

1 

UNCLASSIFIED 



UNCLASSIFIED 


— Aerospace and Electronic Systems- 


The concept of stabilizing satellite attitude solely by gravity 

gradient forces and a passive damper has received considerable attention 

during the past 12 years. Satellites which have employed the method, 

1-4 

however, often have not performed well. Attitude motions of several 
degrees or more are to be expected even under the best conditions. Failure 
of the passive damper always is a possibility. Major, stability problems 
have been experienced even when the damper was functioning normally. 
Particular difficulty has been encountered due to thermal bending of the 
long booms which usually are needed to obtain adequate moments of inertia. 

The yaw motion on triaxially-stabilized satellites has been very trouble- 
some. Many systems which are intended to improve the performance by supple- 
menting the gravity forces by semi-passive, semi-active, or active control 
have been proposed, and some have been built and tested in orbit'* . Reaction 
jets are one of the methods by which performance can be augmented. Emergency 
control of satellites which are strictly passive during normal operation is 
one potential application; this concept should be potentially attractive 
on satellites, such as RAE-B, where a jet system already is required for 
other purposes. 

If active control forces are to be applied to satellites with highly 
nonrigid parts, such as inertia booms, careful consideration must be given 
to the excitation of the elastic vibration modes. At present, there is 
considerable interest in the dynamics, stability, and control of nonrigid 
satellites . The present study emphasizes nonrigidity effects. It is 
applicable to situations in which the purpose of the active control is to 
attenuate the attitude motions without generating unacceptably large 
structural vibrations and also to cases where attenuation of both the 
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initial attitude and structural motions is desired. The work can be applied 
to rigid satellites by deleting the nonrigidity portions of the mathematical 
results. The study is devoted to the basic dynamics aspects of the control 
problem. The related problem of state estimation is given only secondary 
attention. 

Ihe attitude control of rigid satellites by jets has been studied 

in numerous previous papers. A partial listing is given as Refs. 9 to 13. 

The present study differs from earlier ones mainly in that (1) it includes 

nonrigidity of the satellite and (2) it uses the Hamiltonian H to establish 

the control laws. Previous use of H in satellite studies has been directed 

mainly toward its application as a Lyapunov function. The first such study 

was performed by Pringle 1 ^. Meirovitch, Likins, and Budynas and Poli also 

have published papers in this area. H al3o was used by Mackinson and 

Bainum^'* in analyzing the performance of a magnetically damped, rigid 

gravity gradient satellite and by Bowers and Williams^ to optimize the 

timing criteria for the first boom deployment on RAE-I. Although he did not 

17 

use the word "Hamiltonian it also was employed by Watson in a study of 
the capture of a rigid satellite by the gravity field. 

The’ Hamiltonian used here is the H of the nonrigid satellite 1 s 
mechanical state X relative to the rotating local vertical-orbit pole 
coordinate frame R under the approximation that the orbit is Keplerian. 

It will be assumed, as necessary, that the active control forces u are 
balanced so that they do not alter the orbit. Consider first the idealized 
case of a circ ular orbit and negligible disturbances and internal damping. 
Assume that all three moments of inertia are unequal, that the satellite 
includes no parts —such as free rotors— whose motion is not constrained by 
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stiffness forces, and that it is passive except for the control jets. 

When U = 0, the satellite then will possess time-invariant stable equili- 
brium states JC 5e . H can be defined in such a manner that it is zero when 
H then will be positive for any X it also will be time- 

invariant when H — 0 regardless of theX motions. H thus can serve as a 
metric of the displacement of X from X Se . The purpose of the on-board 
control system can be regarded as being to drive H to the smallest possible 
value. (H can be driven to zero only if the mathematical model of the 
dynamics is completely controllable by the thruster systems). H or its 
reduction AH, therefore, can serve as performance measures for developing 
the control laws. The H criteria must be supplemented by a constraint on 
the DC Se about which capture is achieved in applications where the satellite 
is tumbling at the start of the control operation and not all of the stable 
equilibrium orientations are suitable for the mission. Under some conditions 
- particularly a highly nonrigid satellite, large control forces, and/or 
large initial attitude rates - control laws based solely on H can yield- 
unacceptably large structural vibrations; in these cases, structural vibra- 
tion mode amplitudes impose additional constraints on the control. 

H will include cyclic coordinates if two moments of inertia are 
identical or if the satellite includes elements whose motion is not con- 
strained by stiffness forces. Control criteria based solely on H can 
control the velocities of these cyclic coordinates, but not their magnitudes. 
The H viewpoint of active control still is applicable when disturbances, 
orbit eccentricity , or internal damping are not negligible. H then, 
however, will not be time-invariant when u = 0, Disturbances and C can 
produce short term variations in H and, under some conditions, secular growth. 
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Reaction jet systems, particularly those that are operated with low duty 
cycles, should be effective mainly against the secular component . Internal 
damping is beneficial, 3ince it yields H ^ 0. In some situations, active 
control operations yield the side benefit of increasing the fraction of H 
which is contained in structural vibration modes where internal damping is 
present. 

In the following sections, 3x1 Cartesian vectors are indicated 
by an over-bar. Lower case letters with an under-bar are column matrices. 
Upper case letters with an under-bar are rectangular or square matrices. 

A prime indicates a matrix transpose or a row vector. A comma signifies 
differentiation. The double subscript summation convention is employed in 
portions of the work, but only where specifically indicated. Subscripts 
i and j span the range 1 to m where m is the total number of generalized 
coordinates. ^ ^ j and y span 1 to h where n = m-3 is the number of 
"structural" coordinates. P and 6 span 1 to 3 to indicate rotation 
coordinates. spans 1 to 2 m to indicate state variables. o< and /3 span 
1 to s where s is the number of on-board jets. The Nomenclature Section 
includes only symbols which are not defined elsewhere in the paper. Specific 
definitions are not given for symbols which are believed to be self -defining 
from their usage in the equations. 
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BASIC EQUATIONS AND CONTROL POLICY 


This section develops an approach to the problem of establishing 
control laws for active augmentation of the dynamics of nonrigid gravity 
gradient satellites and derives the necessary analytical model. The follow- 
ing two sections employ this material in studies in which the controller forces 
first are modeled as impulses and then as rectangular pulses. The Hamiltonian 

rate H under the action of control forces and disturbances is significant to 

18 ’ 

the work. The classical equation H = H, t is not applicable because it 
does not include the effects of nonconservative forces and because it 
assumes that H is expressed as a function of generalized coordinates q, time 
t, and generalized momentum rather than q, t, and velocity which will be used 
here. To obtain H, consider a mechanical system whose configuration is 
specified by an mxl coordinate vector q. Using subscript summation convention, 
H and the Lagrangian L can be expressed as follows: 


h = -s' <§. 

L - .5M; ( ' U <2> 

where U ~ V — T 0 + K 0 (3) 


The operations required to derive the desired H equation are (a) substitute 

18 

Eq. (2) into Lagrange's equation , (b) multiply this result by and sum, 
(c) form H by differentiating Eq. (1) and simplify using the result of step 
(b), and (d) convert to matrix notation. The result is: 

h = j-'lf? - f 3 + <4) 

where is the mxl generalized force vector. 
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In the present satellite application, £ defines the system's con- 
figuration relative to a rotating local vertical reference frame R. The 
origin of R is the satellite's c.m. . The orbit will be constrained in the 
modeling to be Keplerian. Equations for K, b, and U are given in Appendix A. 
The partial time derivative terms in Eq. (4) arise from (1) active on-board 
alteration or control of the satellite's geometry or mass distribution and 
(2) disturbance phenomena that are not included in fq. Often, it is easier 
to model actively-controlled geometric variables as independently-derived 
time functions J^(t) than it is to include them in a- In this case, the 
partial t derivative terms in Eq. (4) can be the main mechanism through which 
active control of H is accomplished. Such systems, however, will not be in- 
vestigated in the present study. The work instead will encompass only 

18 

scleronomous-coordinate models in which the active control is generated 

through £q by point forces U ^ and which otherwise are passive. Assuming 

that the satellite's mass properties are affected insignificantly by the 

expenditure of controller fuel, M, t in Eq. (4) then will be zero. U, t and 

b,t, when nonzero, will affect H as disturbances. Eq. (A. 3) shows that 

b, will be zero if orbit eccentricity e is zero. U,. will be zero if e = 0 
— t z 

and if, in addition, phenomena, such as temperature gradient variations, that 
make V 3 vary explicitly with t are negligible. 

fq consists of components fq ds fq^ and f^ due to disturbances, 
internal damping, and the active control forces respectively. In the present 
application, fq dg can include not only non-conservative forces, but also 
conservative phenomena such as gravity gradient harmonics which act as dis- 
turbances and hence should be excluded from H. Let the satellite contain s 
jets and let u (t) be the s x 1 vector comprised of the scaler values of their 
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force outputs. (In the present terminology, a unit which can apply forces 

of either polarity is considered to be a single jet). Using the classical 

18 

equation for generalized forces , and u can be shown to be related 
through an equation f^ = Y u where Y = Y (g ) . Let Qr^ be the s x 1 vector 
of the scaler values of the translational velocities of the jets, relative 
to frame R, along their force axes. It can be shown that and 4 axe re ~ 
lated through Y ’ . 

Ofa = V i (5) 

Lumping all H components due to disturbances and orbit eccentricity into a 

• • • 

single term H^ g and letting and H dm be the components due to active control 
and internal damping respectively, Eq. (L) now becomes 

H = H + H + Hj (6a) 

U rim ds 

where Hy ~ /Vlf ^ H u . - (6b) 

The time integral of H u is the work done, relative to frame R, by u. A 

block diagram of the system is given as Figure 1. 

The purpose of the active control can be regarded as being to reduce 

H and/or to maintain it small. The present paper will employ as a basic 

control principle the tenet that none of the individual jets ever should act 

directly to increase H. When any jet ®< is fired, the correct polarity of 

its force then is opposite to the instantaneous velocity . ihis 

policy is not optimal in all problems. The solutions of minimum time problems, 

11 

such as the rigid gravity gradient satellite one considered by Zach , fre- 
• • 

quently yield H > 0 at times during the operation. The advantage of the 
present 0 approach is that it provides a relatively simple basis for 

developing control laws in nonrigid satellite problems where the dynamics 
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may be nonlinear, high order, and difficult to predict accurately. De- 
pending on the thruster-satellite system and the operational requirements, 
many control policies are possible within the framework of this 0 

criteria. Eq. (6b) shows that the effectiveness of any jet in reducing 
H is directly proportional to the magnitude of the instantaneous 

velocity . In applications where fuel utilization efficiency is the 
main requirement, a simple on-off policy consists of (1) turning each indi- 
vidual jet c( on at times when its \ / VZf c < I is close to a local maximum and 
(2) de-energizing o< when I has been reduced to zero or a specified 

nonzero value. The policy can be modified as necessary to incorporate con- 
straints on the control, such as employing the jets in pairs. The near- 
maximum | hr^l criteria for energizing the jets is abated in applications 
where rapid H reduction is more important than fuel economy. By driving the 
nrijA 's to zero a sufficient number of times it can be anticipated that, in 
the absence of disturbances, a condition will be reached wherein the between- 
pulse values of the /ir^' s remain negligibly small. This signifies that any 
remaining motions of the satellite are in uncontrollable modes, and no farther 
reduction of H by jet firings is possible. 

The studies to be presented in the following two sections require 
dynamical equations for the satellite's response in q. , g_ , and Qr^. The 
generalized coordinate method which will be used was employed previously in 
Ref. 19. The m generalized coordinates are of two types, q^ to q^ are a set 
of Euler angles 9 which define the orientation, relative to frame R, of an 
arbitrarily-selected body frame b of the satellite. The remaining n = m-3 

coordinates specify the satellite's configuration relative to frame b. 

.... 20 

The term "hybrid coordinates" has been coined by Likins and Wirsching to 
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indicate "mixed" coordinate sets such as the L Q\ r^J one. The portion of 

£ is omitted if the satellite is assumed to be rigid. In the equations of 

motion, it will be convenient to use the angular velocity vector CO ^ of frame 

• • 

b relative to frame R instead of © . and © can be related through a 
transformation CO b =J (jJ @ © where I (u6 ,=I oe (©) . Let arbe the m x 1 velocity 
vector nr = lo^ / r^J. nr and q then are related through nr- z ^ where 


T - 

lf)r f 


~J^OOQ .O 


£ 1 

Several rectangular matrices in the dynamics equations must be transformed 
when switching between nr and q velocities. Matrices that are employed with 
nr- will be indicated by a subscript nr . In QT velocities, Eq. (5) becomes 




( 7 ) 


where 


- 1 


Y' - y f t 

/ ' r_ J. — V ~C^ 


J.nr 


The kinetics model which will be used was obtained by modifying 
the equations developed in Ref. 19 to make them directly applicable to the 
present application. The resulting equation is 

hA^ 0j~ N't, r ~ U -f- -i P 2 (^) 


Appendix B presents detailed expressions for the terms appearing in Eq. (6). 
Similar-appearing equations were used in flexible satellite studies reported 
in Refs. 20 to 22. The equations in these studies, however, all were lineari- 
zed, limited to satellites comprised of a rigid central body to which non- 
rigid au xilli ary bodies are attached, and also limited to very specific 
representations of the nonrigidity. These restrictions have n ot been imposed 
on Eq. (8) and the associated equations in Appendix B. The advantage of the 
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use of £T rather than g velocities is largely due to the fact that the 
method reduces the dependency of the rectangular matrices in Eq. (8) on Q. 

and Y^ are functions of f]^ . N v is a function of and ©. Unless the 
satellite is highly nonrigid, however, it usually should be possible to omit 
the l\ dependency in Y^, and . Mat is symmetric and positive de- 

finite. N^ar results from (l) gyroscopic forces due to and (2) viscous 
damping. is skew symmetric if damping is negligible, and £ 2 include 

both actual forces and also "apparent forces" due to kinematic effects, 
is comprised of conservative forces due to internal stiffness and central 
force field gravity gradient; it also includes the centrifugal potential force. 
f 2 consists of apparent forces involving quadratic terms in nr-j it also in- 
cludes disturbance forces. 

IMPULSE MODEL OF CONTROL FORCES 

This section will use the material developed in the preceding section 
to investigate the case where the control forces are of sufficiently large 
magnitude and short time duration that useful information can be obtained by- 
modeling them as impulses. Ihe advantages of the impulse representation are 
that it is the easiest model to investigate and it yields the simplest 
result. The study does not require linearization of the equations of motion. 
The weaknesses of the impulse assumption will be discussed at the end of the 
section. Let subscripts 0 and 1 indicate conditions at the start and end of 
a control action which consists of the simultaneous firing of one or more of 
the s on-board jets. Let Uj be the s x 1 impulse vector. Uj is the time ; 
integral of u. Uj will produce step changes A arand AH . Since g undergoes 
no change during the impulse, the matrices M^, , and Y^ in Eq. (8) will 
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be constant. f^, and f^ do not affect the A at step, 

can be shown to yield 



Zl nr 

— Mfir' Y [A 


AH - 

Tuo +* 

where 

D ~ 

v; Mr Y«- 


5 * Uf 0 & j c 


Eqs. (6) and (8) 

(9a) 

(9b) 

(9c) 


The step change in^r u is = Dup Any element of D thus is the velocity 
step change along the force axis of jet which is generated by a unit im- 
pulse from jet /3 . It can be demonstrated that D is symmetric and at least 
positive semidefinite. 

& 

Eq. (9b) now will be used to derive the impulse vector u£ which 
yields the minimum (i.e., the most negative) value AH* of A H. It will be 
assumed that no amplitude limits are placed on any of the u l£a< 's. The 
derivation of u^* then constitutes an unconstrained static minimization 
problem 2 ^. The solution given in the following paragraph will demonstrate 
that Uj* is unique if and only if D is nonsingular. D will be only positive 
semidefinite, and hence singular, if the satellite possesses redundant jets. 
The analytical difficulty which this causes can be overcome by imposing a 

selected set of. 3d constraints on the relations between the u^ 's where s^ 

25 c c 

is the degeneracy of D. A constraint equation u j = C , therefore, 

will be assumed. Uj C has dimensions s c x 1 where s c — s— s^. The constraints 

might, for example, be set up to restrict the 3x1 Cartesian vector of the 

total impulse to be zero; this eliminates the potential problem of altering 

the orbit significantly by the jet firing. The s c x 1 impulse vector u I ° 

can be regarded as the output of a ficticious set of s c jets which will be 

c* 

referred to here as the "constrained jets." u^. will be unique, since the 
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I c 

transformed s c x s c matrix D c = C c pC c which appears in Eq. (9b) when u^. is 

employed is nonsingular. For the present purposes, it will be convenient to 

impose an additional set of constraints on the jet firing. The discussion 

in the preceding section indicates that use, in any single firing operation, 

of all the s_ constrained jets often should be undesirable because it utilizes 

V 

the available fuel inefficiently. Let u a be the s a x 1 impulse vector from 

8 / 

the s a constrained jets which actually are used in a given firing, u^. and 

c , ^ a 

Uj can be related by an equation u ^ = b U 


.ca 


has dimensions x s 
c a. 


where s a < s c . Each row of C ca will contain either (1) all. zeros or else 
(2) all zeros except for a single element with a value of unity. 

Substitution of the relations specified in the above paragraph into 
Eq. (9b) yields 


~ Ofu o W j 


6 & Pa 


(io) 


where /ir' a = C 1 r\r , D d = C 1 D C, and G = C c C ca . u“ now can be 

obtained by differentiating Eq. (10) with respect to u a , setting the left 
side to zero, and solving. It also is necessary, in principle, to differ- 
entiate with respect to u a a second time in order to check that the resulting 
A H* is a minimum rather than a maximum. The result is: 

.a 


a* - 1 

u_ — -D /ir~‘ 
—I -a 




p FI * — 3, 

a h *= .5 orjo ^ 


d 2 AH = 
.2 a 
d % 




uo * 


D 

-a 


(11a) 

(lib) 

(11c) 


Since is positive definite, AH* is a minimum, 


24 
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u* can be obtained from u^* using u^ = C u^. . The equations 

dP' 2 

given in this section can be used to show that yields ^ — * Thus, 

a *VC 

the optimal impulse u^. is the one which drives the force axis velocities, 
relative to frame R, of the s constrained jets to zero. Eq. (lib) shows 
that the greatest reduction in H in any single firing operation is obtained 
by using all s c constrained jets except any whose initial velocities /Tr^ 
are zero. Eq. (lib), however, also shows that the AH* component which is 
contributed by the optimal impulse u a "^of any individual constrained jet ^ 
is proportional to the magnitude I ^ I of its initial velocity. Thus, from 
the standpoint of fuel efficiency, the most effective firing procedure is to 
fire each constrained jet only at times when its/'V^Jis passing through a 
maximum. 


a* 


Practical satellites usually should be sufficiently rigid that the 

variation of IP*' with t\ can be omitted when establishing the criteria for 

a * 

active control. Eq. (11a) shows that u^ then is not dependent on the 

A 

complete state X of the satellite. Rather, it is dependent only on Q£ uo 
If the satellite's sensor system is such that a measure of nr^ can be obtained 
without generating an estimate of X- , actual knowledge of _X then is not 
needed for control. 

While the unlimited-amplitude impulse model of the jet outputs 

yields a simple criteria (Eq. (11a)) for determining the instantaneous optimum 

control effort, the technique encounters two major difficulties when the 

satellite is highly nonrigid. First, even if the mathematical model of the 

3 * 

satellite and jet outputs were perfect, it is still possible that the U x 
computed by Eq. (lla) might generate unacceptably large structural vibrations. 
Placing amplitude limits on the impulses is one approach toward alleviating the 
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potential problem. With this method U T would be obtained by solving 
Eq. (10) subject to control constraints / ^ ^ x A wax ^ ^ ~ ^ ° c ^ ' 

A less fallible, but more difficult, approach is to place the amplitude con- 
straints directly on the magnitudes of the structural vibration modes. This 
method is used in the following section in the rectangular pulse approximation 


study. 

The second difficulty concerns the compatibility between the approxi- 
mation that the jet outputs are impulses and the implicit assumption that, 
for the present purposes, the nonrigidity of the satellite can be modeled 
adequately by a finite number of coordinates y . The use of a finite number 
of (\ v ' s to model a continuous structure requires, in principle, the truncation 
of an infinite set. The ^,,/g which are included are ones which contribute 
mainly to the lower frequency modes of vibration; theft's which are omitted 
contribute mainly to the higher modes. Before Eq. (11a) could be used in an 
actual nonrigid satellite application, it would be necessary to verify, by 
supplementary analysis, that the excitation of the neglected vibratory modes 
by the short force pulses will not be great enough to cause difficulties. 

Also, Eq. (11a) will be accurate only if the time durations of the actual 
pulses are short enough in relation to the periods of the vibratory modes 
which are included in the model that the impulse approximation is acceptable. 
The filtering action which tends to attenuate the effect of non-impulsive 
inputs on higher vibration modes is not encountered in the present model. The 
computed therefore, depends on the number of rj y 's that are 

used. 
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RECTANGULAR PULSE MODEL OF CONTROL FORCES 


In the present section, the active control forces will be approximated 
as rectangular pulses. This representation alleviates, to a considerable 
extent, the impulse model difficulties noted at the end of the previous 
section. The on-off control policy discussed earlier will be used. The 
problem which will be considered is the computation of the optimum pulse 
duration^ when the "height" of the force pulse and the satellite's dynamic 
state Iti at the start of the pulse are known. In an actual application, the 
pulse durations would have to be determined a priori by such a computation if 
lags or other deficiencies in the sensor system prohibit a "closed loop" 
determination of the proper jet shut-off points using real-time information. 

An analytical investigation using non-impulsive control forces re- 
quires, as a practical necessity, a linearized, constant coefficient model of 
the dynamics. The first step, therefore, will be to linearize Eq. (8). A 
model which is applicable at least over the time spans of the individual pulses 
is desired. Linearization of Eq. (8) limits the study to cases where the 
quadratic nr terms in are negligible over the duration of a pulse. The pre- 
sent work will omit disturbances. then can be deleted entirely. The study 

also must be restricted to cases where the coordinate changes Ei o, during a 
control pulse are sufficiently small that M ^ and can be considered 

constant and that (~ 0 can be approximated by the first two terms in a Taylor 
series expansion. Most satellites should be sufficiently rigid that no major 
restrictions on the model's validity are imposed by omission of the higher 
order effects due to the dynamics of the structural coordinates H . However, 
when the satellite's initial attitude is not close to the stable equilibrium 
attitude, the subsequent Q motions can be large enough to severely degrade the 
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linearized model's validity. In such cases, the linearized model is applicable 
only if the control pulses are sufficiently short that only small changes in 
£ are experienced during their . dur at ion. The study then must be limited to 
large amplitude-short duration pulses such as those which were modeled as 
impulses in the preceding section. 

In linearizing Eq. (8), will be transformed into a velocity vector 
which can be approximated as the time derivative of an m * 1 coordinate vector 
$ . Consider a structural reference condition P( which is close to the 
initial value £] • Let $ 4 . to be chosen as p. The condition 
encountered when the satellite is in equilibrium usually should be a suitable 
choice for . In establishing the rotation coordinates to , use 
will be made of an auxilliary reference frame "r " which has zero angular 
velocity relative to frame R and which, at t , has an orientation close to 
that of frame b. Let & . to ^ be a set of Euler angles which specify the 

attitude of b relative to r . ^ = 0 occurs when b and r are aligned. 

• — ... • - . 

and can be related through an equation G ^ ff where 0 - 7^ ( 

• ~ 0 

and (o') Similar ly, Ar and jS can be related through nr -= JT^ @ where 

7 ^ and ( 0 ) - J ; T rr& is defined analogously to used 

previously. Employing this relation in Eq. (8), premultiplying by_T^ and 

9 0 

dropping r £ and a quadratic velocity term involving 7^ 6 yields 


Jyir# — ~LrtT0 + J/tr# jfw M 


Eq. (12) will be linearized about the reference ( r ) condition 
= 0. Superscript r will be used to indicate quantities which are evaluated 
at 0 = 0. When expanded about /£ = 0 and linearized, the left side of Eq. 
(12) reduces to (i_ + j\]^. Q . It can be shown that P D a 

18 


UNCLASSIFIED 



UNCLASSIFIED 


— Aerospace and Electronic Systems 
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where the symbol LU indicates L 0^ J • Expanding in a series 

about /2 = 0 and retaining only the first two terms yields U, ^ - 
■ ; 

The elements of are listed in Appendix B. IheJ^j premultiplication in 

j. r~ 

Eq. (12) was performed so that the resulting stiffness matrix K^would be 
symmetric. Expanding Xir /3 Y»r d and retaining the only the first two terms yields 
u — P 1 "^ where the columns of the unsymmetric m matrix P^are of 
the form d * Det ai led expressions for the jO£ s are given in Appendix 

B. In applications where U is constant over the time intervals that the 
linearized, constant coefficient model is to be used, P can be included by 
adding it to k^- as an additional stiffness effect. It is not certain that 
P_ h will be negligibly small with respect to kjlrin all cases. The first three 
P h rows appear to be potentially comparable in size to the similar k/v- rows 
whose elements are nonzero solely because of gravity gradient and centrifugal 
forces. It is believed, however, that omission of P h 0 usually should be 
justifiable largely on the grounds that is much smaller than U . Conse- 
quently, Jf V will not be included in the remainder of the paper. 

Making the above-noted modifications to Eq. (12) and performing 
similar operations on Eq. (7) yields the desired linear equations. 

3 = fd i (13a) 




f 1 'JV.'l 


In order to solve Eq. (13b), conversion to state variables X is advantageous. 
Letting X - L& \ (b J ^Eq. (13) can be converted to the form 

^ - £ C ' I °) 1 (Ua) 
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X - A a + § 


- I 


<- 


1 








o 


. i 





o 


(14b) 

(14c) 


(14d) 


Eqs. (14) next will be transformed into modal coordinates This is con- 
sidered desirable in the present work because it enables a simpler and more 
informative solution to be obtained. The 2. m eigenvalues > ^ of A will be 

assumed to be distinct. Let X“ ^ inhere the am columns A^ of S are the 

r -i — 

eigenvectors of A. Let R be S , and let K indicate the columns of R. Each 

5 k' ; ' f 

will be separated into two my/ components j. this same 

notational technique will be used later with g and , 
will be indicated by ft* I dimension vectors rz£ 

(14b) then can be written in the form 


The s columns of 
, Eqs. (14a) and 


Z. tyj 








. jj 


V- 


m -- 


(15a) 


(15b) 


The solution of Eqs. (15b) is 

r J 


§ 


- <z 


(16) 

"Ci 
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More than one jet may be in use at any given time because of their 
non— zero firing durations and because of possible constraints on their oper- 
ation - such as firing them in pairs. Subject to restrictions which such 
constraints may impose, each jet should be energized when its j 
becomes appropriately large, (The present terminology assumes that each jet 
can apply forces of either polarity). o<( is de-energized at, or slightly 
before, the time that changes sign. Assuming that the jet pulses are 

rectangular, g in Eq. (16) will be constant during "constant control con- 
dition" ( CCC) time segments during which no jets are turned on or off. 
Inserting the constant g condition into Eq. (16), the g responses during 
a CCC are determined to be 


Jrs C C -L°*f 8 )z 


Vt-t. 


N 

OJ 


\ n 

- 


(17) 


Assuming that the CCC is ended by de-energizing a jet because its terminal 
f\rr , condition is reached rather than by energizing a new one, the CCC end 

14 & 

time can be computed using Eqs. (15a) and (17). Ihe computation consists of 
a determination (by iteration for example) of the first time point at which 
the cut off f\f , ^ is reached by any of the jets that are in use. An equation 
for the due to jet action during the CCC can be derived using Eqs. (6), 
(15a), and (17). Ihe result is 


i:'± i ns) 




AH' <4 




r-y 




f 
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As was noted previously, a drawback of the pure H reduction criteria 
for establishing control laws is that it potentially can generate unacceptably 
large structural vibrations even if the mathematical model can be considered 
perfect. On some satellites, this problem can be alleviated by employing the 
thrusters in sets so that the sensitivity of certain of the modes to the 
firing is reduced or eliminated. The present paper will consider the supple- 
mentary technique in which the jets are turned off prior to reaching the 
^ = 0 condition when necessary to prevent the magnitudes of critical 
structural modes from exeeding selected limits. This approach encounters 
some difficulty if the satellite is highly nonrigid and the attitude motions 
following the control application are sufficiently large that the subsequent 
structural dynamics are significantly influenced by nonlinear effects. The 
problem is that energy transfer due to nonlinear coupling, in principle at 
least, can cause the amplitudes of some of the vibratory modes to increase 
after the jets are shut off. This phenomena, however, will not be considered 
in the present work. The study, instead, must be restricted to cases where 
the previously— developed linearized model is adequate. 

It will be assumed that, for the present purposes, internal damping 
can be considered negligible and that the 2n eigenvalues ^ of the structural 
modes thus will be imaginary: \ f will be complex. 

Since the modes occur in complex conjugate pairs, only n need be included. 
Using the rectangular control pulse approximation, Eq. (15b) can be integrated 
once to yield the following equations for the real (R) components of the 
modal motions. 
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f ( liVtZ ~ ‘ <V (19a) 

where ^ V R “ ^ VP ^ ^ a V X ' p (19b) 


^yu ~ 

-± 

Jfavx'gau 

S 

\ 

(19c) 

/\s 

i 

7)vr. 


(19d) 


Kyis a positive constant which depends on the initial conditions. o 3cf and 

are the components of due to U and ^ . Jfayjl is the imaginary com- 
ponent of Jfav ♦ Ihe trajectories defined by Eqs. (19) are circles in the 

/V ~ . /\y 

oTy^- phase planes. When M = 0, the centers of the loci are at 

= 0. Firing one or more jets shifts the centers to ^y R ~ — <> 

A suitable choice for the amplitude limit constraint equations is 


~ ~ 2. / , 2. 
^ Vs ^ L >' 


( 20 ) 


It is assumed that an modes V are inside their boundaries at the start 

of the control operation. Figure 2 shows a typical response. If, as the 
result of a control action begun at t Q , a mode V reaches its L~. y boundary at 
t^, a modification to the control is made at t^ to prevent l~^> from being 
crossed. This normally would consist of de-energizing the jets that excite 
mode V . It would take precedence over the 0 criteria. The point 1 

where V reaches l-y can be determined by solving Eqs. (19a) and (20). The 
result is 



( Ly + ^ yu - Ky ) / 2. yu 


(21a) 


/V 

'IT* 


VRl 


~ ± ( L‘ 


7t VR. I 


) 


(21b) 
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An imaginary ctr yR( indicates that the L y boundary will not be reached. The 

time t, - t can be shown to be 
1 o 


(VCO - 


2 - 


- f 

Sin l 


>b_ 


r r a/ < 

*“* VIZI ~ l)v/Zo ) 


( 22 ) 


+ ' Ar vtzoJ J 


. ? 


) 


CONCLUDING REMARKS 

The Hamiltonian provides a useful measure of the displacement of 
gravity gradient satellites from equilibrium. It can be used as a tool for 
developing supplementary active control laws for rigid and nonrigid satellites. 
Such approaches deal with the composite system in a unified manner and can be 
used in cases the more common methods are difficult to apply. The main draw- 
back appears to be the potentiality of exciting intolerably large structural 
vibrations when the satellite is nonrigid and the initial attitude motions 
are large. The velocities A^of the thrusters along their force axes are of 
primary concern for control. The impulse model of the jet outputs yields a 
simple and practical control criteria when the forces are large and the sate- 
llite is near-rigid. This approach must be applied with considerable caution, 
however, if the satellite is highly nonrigid. The rectangular pulse model of 
the jet outputs has far wider applicability, but leads to more complicated 
control law mathematics when an "open loop" computation of the jet cut-off 
times is needed. 6 
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Figure 2. Phase Plane Response of Mode y 
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APPENDIX A 

EQUATIONS FOR THE HAMILTONIAN OF NONRIGID GRAVITY STABILIZED SATELLITES 


Ihis Appendix and the following one are included to aid in inter- 
preting and employing the material in the main body of the paper. Ihe equa- 
tions for H can be derived straightforwardly using basic principles. The 
general form of H was given in Eqs. (l) and (3). Eq. (1) can be converted 
into matrix notation and can be expressed ir. Or velocities by the transform- 
ation /tr"-T \r<^ %- • In the present application, the terms in Eqs. (1) and (3) 
can be shown to be 



r jy 

-A"' 

z t 

1 

~A 

A 1 

— 



(A. la) 


M- 

U - 


-T- f ~T 

1 rrc/ fir * fir a, 

- V - ~ r (A. lb 

- 2yii mn e cos Y -- 

a-Ci-e 2 ) 3, £3e K J e R -Ti-jdc J1( A2 ) 

-- 6 W s 5 } 'a'j, l ■+ Zb' A} l T K u 


The term b which is used in Eqs. (2) and (4) is 

b : - cs/cs i- a 1 ) 


k - 


'fir 


A-l 


(A. 3a) 
(A. 3b) 
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and A^ arenxnand 3 x n respectively. IVie elements of Ap" and columns 
j, 2 of A 2 are obtained by integrating functions of e. and J 7 over the composite 
satellite; the exact definitions of A^,y and A^ are given in Ref. 19. Jl is 
comprised of coordinates which are treated as independently-derived time 

' 0 

functions and hence not included in When the satellite is passive, 

2 — 12 

is zero. A - ^ and A ^ are defined similarly to A and A . 


A-2 
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APPENDIX B 

EQUATIONS OF MOTION OF NONRIGID GRAVITY STABILIZED SATELLITES 


This Appendix lists detailed equations for the terms in Eqs. (8) and 
(13). Eq. (8) was obtained by converting the basic nonrigid satellite equa- 
tions in Ref. 19 into a form which is more convenient for the present study. 

The main modification was the use of angular rate C relative to frame R in 
place of to relative to a nonrotating frame. to 3 to^ i and their first time 
derivatives with respect to frame b are related by: l o —to^~t : tv ~ ( v^ ■+ <(<3^ e 

The other modifications are mainly of a notational and algebraic manipulation 
nature. It might be noted in particular that the symbol is used here in 
place of Ref. 19 's q. The main restrictions on the equations given in this 

O 

Appendix are that C and^f are assumed zero. 

/VI ^ was defined in Appendix A. The other terms in Eq. (8) are 




(B.la) 


/\J f/ 2 J - 1 T hace J ^ • L Vfc 3 


(B.lb) 


K ty ■=■ -2 Ay* “>R 

+ N lluy (B.ld) 


\/ 

I f\j~ 



(B.2a) 
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W,/rf - ' Z UA 



& C •“ 3 6o| O' J “’r 


(B.2b) 


(B.2c) 


(B-3a) 


(B.3b) 


r 5 ^2. ? (B-.4a) 

_ r_ W 

^ = S ds - hy ' & 6 2 H + ^ ^ (B - 4b) 



$ ““ 3 ^ 


|\J y and Zy are the columns of and Z. The double subscript summation 
convention is used with and V . Superscript S indicates that the Cartesian 

vector is arranged in its 3^3 skew— symmetric form. The A terms are defined 
in Ref. 19. is an n x 1 vector formed by differentiating U with respect to 

the elements. Vectors which appear in the first three rows of Eq. (8) are 
resolved on frame b. 

The elements of K^in Eq. (13) will be derived next. Limiting the 

* 

work to the <Z - O 9 - _£ case, Eq. (A. 2) can be written as 

U - V s ^fz ^ 'H. ^-bf- <£r 

ZZ b * r f -r-b — h -r- b •. (»5) 

" ^f> irb h _ - I hdce i 3 f l 


Superscripts b and r indicate the resolution frame. C.^)- is the direction 
cosine matrix from \r to b. is constant, C-bb' = '^-bi- — J f H 3 5 

B-2 
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and V c - V s ( n ^ • since /£ — I?' I J i we thus have ( , 

As was noted earlier, «^ - [_ .Separating l<^ into submatrices 

yields 



^ a> (p 

^r><£A 




(B.6) 


The elements l^y of the n x n matrix U, ^ can be obtained directly by 
differentiating Eq. (B.5). The columns of the 3 x n matrix ^ can be 

P f 

shown to be equal to - G- 0jV • will be unity, since it is to be 

evaluated at the reference condition (p = 0. G~ 0 y is formed by differentiating 

Eq. (B.3b) with respect to f) y j J" is the only term in Eq. (B.3b) which is 
a function of fl . The elements /frof the 3x3 matrix are obtained by 

differentiating Eq. (B.5) with respect to <Pj and <Pp . As an aid in the 
derivation, it is advantageous to first specify the sequence of the 0 rotations. 
Let the order in going from frame r to frame b be 0-^, 0£, 0^ where the sub- 
scripts 1, 2, 3 also indicate axes X^X. 1 ^ , It then can be shown that 

cJ bv 6P - P6 - -C 6 ; where ^ is the larger of the 6 pair 

and /s? is the smaller; the C^'s are unit vectors along the instantaneous 

axes of Euler angle rotation. Using this relation and Eq. (B.5)j ^5 p(> can 

be shown to be 

v ^6 - - ( 4 K{i C3 e* ( 2 + U (B . 7) 
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generated when expa: 

h 


P )r -r \J tJ 

generated when expanding _[nr/i J. v ' — - 
in Eq. (12) will be considered next. Its columns ^6^ are 


A L fg 
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where & u £ is the torque component along axis . The remaining ^’s.are 

V 

U 

(B.10) 
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